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Abstract 

Bertrand's theorem in classical mechanics of the central force fields attracts us 
because of its predictive power. It categorically proves that there can only be two 
types of forces which can produce stable, circular orbits. In the present article 
an attempt has been made to generalize Bertrand's theorem to the central force 
problem of relativistic systems. The stability criterion for potentials which can 
produce stable, circular orbits in the relativistic central force problem has been 
deduced and a general solution of it is presented in the article. It is seen that 
the inverse square law passes the relativistic test but the kind of force required for 
simple harmonic motion does not. Special relativistic effects do not allow stable, 
circular orbits in presence of a force which is proportional to the negative of the 
displacement of the particle from the potential center. 

1 Introduction 

The central force problem in non-relativistic classical mechanics is one of the most useful 
topics in physics. Closely linked with the central force problem is the Keplerian orbit 
theory which is a cornerstone for understanding planetary motions in the solar system or 
motion of electrons near the nucleus. In classical mechanics there is an important theorem 
called the Bertrand's theorem which proposes that there can only be two types of central 
potentials, the Coulomb type and the simple harmonic type, which can produce stable, 
circular orbits for particles moving around the potential source. A good presentation of 
the Bertrand's theorem can be found in Ref. [Tj. The present article tries to generalize the 
results of Bertrand's theorem when the orbiting particle can have relativistic velocities. 

In this article we first set up the relativistic orbit equation for a particle in a central 
potential presumed to be dependent on the radial coordinate only. The relativistic central 
force orbits were previously studied in Refs. (2J HJ |5] . A brief description of the central 
force problem in a relativistic setting in a Coulomb potential was presented in the book on 
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classical theory of fields by Landau and Lifshitz [6]. Before one starts the main analysis 
about the stability of orbits of relativistic particles in a central force potential it is better 
to specify the assumptions one makes in arriving at definite results. In the present article 
we use the same assumptions and the approximations as utilized by Boyer in Ref. [2] and 
Landau in Ref. [5]. In the specific references cited above, none of them present a Lorentz 
covariant treatment of the relativistic central force problem. The main reason being that 
all of them assumes a central potential V(r) where r = |r| is the distance between the 
source and the orbiting particle. The form of the potential only depends on the position 
coordinates of the orbiting particle. The form of V(r) is not Lorentz covariant. In such 
cases the results of the whole analysis is valid in a particular frame where the origin of 
the coordinate system coincides with the potential center. 

The references cited above assumes the particle which produces the potential V(r) to 
be static in the specific coordinate system utilized by the observer. If the source of the 
potential does not have any velocity then the retarded nature of the interactions, owing to 
the finite velocity of light, does not complicate the calculation of the orbit of the relativistic 
particle. A specific example will make the point clear. In classical electrodynamics if the 
source of the Coulomb potential V(r) moves with a velocity v s and the orbiting particle 
has a velocity v then the potential V(r) gets a relativistic correction. The magnitude of 
the lowest order relativistic correction to the Coulomb potential was calculated by Darwin 
in 1920 and it looks like 



V{r) 



r)(v 



For a better understanding of the Darwin correction one can look at Ref. |7j|. In our case 
v s = and consequently there will be no relativistic modification of V(r). More over 
we do not consider any general relativistic effects due to V(r) into account. We briefly 
comment on the general relativistic generalization of the central force problem in section 
HI In the present article the background space-time is assumed to be flat. 

In the article it will be shown that the stability condition of the perturbed orbits 
around a stable circular orbit gives rise to a non-linear differential equation for the central 
potential. The Newtonian or the Coulomb potential satisfies the resulting differential 
equation with some restrictions on the possible value of the angular momentum of the 
orbiting particle. Except the Newtonian potential solution we present a more general 
solution of the differential equation for the potential which can give rise to stable, circular 
orbit for relativistic particles. This solution gives rise to a force which is not common 
in physics except its Newtonian inverse square law limit. The equation of the orbit of a 
relativistic particle in such a non-trivial force shows that the orbit will precess and the 
precession angle can be calculated. 

Unlike the non-relativistic case, in the relativistic case there exist no radial effective 
potential minimizing which we can obtain the radius of a circular orbit. In the relativistic 
case a first order perturbation from a circular orbit is enough to determine the stability 
criterion of the orbit. In the non-relativistic case one uses higher order perturbations from 
a circular orbit to specify the form of the potential. In the relativistic case the general 
solution of the form of the potential from first order perturbation from circular orbit is 
such that all higher order corrections becomes irrelevant. As a consequence of this fact 
the general form of the potential which can produce stable, circular orbits for relativistic 
particles contains more parameters than the corresponding expressions of non-relativistic 
potentials. 
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The material in the article is presented in the following manner. The second section 
sets the conventions and derives the orbit equation of a relativistic particle in a central 
orbit. Section [3] generalizes the Bertrand's theorem for the relativistic case. In this section 
the stability condition for the circular orbits will be interpreted as a non-linear differential 
equation for the potential. The solutions of the non-linear stability equation will also be 
derived in section [3j In section H] the connection of the present work with some related 
works which were existing in the literature are discussed. This section gives a wider view 
for the readers who really want to understand the stability of orbits in special relativity 
and general relativity. The last section |5] summarizes the important points presented in 
the article. 



2 The orbit equation 

In this section we derive the orbit equation of the relativistic particle in presence of a 
potential V(r) which is purely a function of the radial coordinate. The Lagrangian of a 
relativistic particle of mass m in presence of a continuous radial potential V(r) is 

C = -mcVl-f 2 /c 2 -y(r), (1) 

where the velocity of the particle v in plane polar coordinates is given as 

v = re r + rOeg , 

where e r , eg are the mutually orthogonal unit vectors along the radial and the angular 
directions. The form of the Lagrangian in Eq. (CQ) immediately shows that the angular 
momentum 

dC 

L= - J = mrV, (2) 
ad 

is a constant, where 

7 



y/l ~ V 2 /c 2 

The total energy E of the particle in presence of the potential V(r) is 

E = mc 2 j + V(r) . (3) 

Although from the definition of 7 it looks like that it is a function of r, r and 9 but it can 
be shown that in a central force field 7 is only a function of the radial coordinate r. The 
reason for such behavior of 7 can be understood from the following reason. As energy 
and angular momentum are constant functions of r, f and 9 we can use the conservation 
conditions of E and L to re-express r and 9 as functions of r, E and L. As E and L are 
constants so in a central force field r and 9 are functions of r alone. Consequently 7 is 
only a function of r. In special relativity the energy of the particle in a central potential 
can also be written as 

E = ^p 2 c 2 + m 2 c 4 + V(r) , (4) 
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where p = |p| 



and 



p = rrryv = m7(fe r + rOeg) 

= p r e r + Pee e , (5) 



p r = rwyr , pg = mjr9 



r 



As because (p r /pe) = (f/r9), we have 



Pr 



L dr 
^dd 



With the above information on the various momentum components we can now rewrite 
Eq. (UK as 

Instead of r we use the variable u = - in terms of which Eq. (E} becomes 

(E - V f = L 2 c 2 V + u 2 L 2 c 2 + m 2 c 4 . 
\dv J 

If we differentiate the last equation with respect to 9 and then divide the resulting equation 
by du/dd we obtain the desired equation of the orbit of a particle of mass m possessing 
momentum p moving in the presence of a general central potential V(r) as 

(V-E) dV 

d9 2 L 2 c 2 du' [) 

Using Eq. ([3]) we can rewrite the above equation in the form 

d 2 u mj dV , . 

de 2 - + u = 'l 2 ^- (8) 

Writing L = j£, where i = mr 2 6 is the non-relativistic angular momentum, the above 
equation in the non-relativistic limit (7 — > 1) transforms exactly to the form we get in a 
conventional non-relativistic treatment of the problem as given in Ref. [T]. 



3 Circular, Stable closed orbits 

Lets define 

Suppose Eq. ([7j) admits a circular orbit of radius r = l/u . For small perturbations 
around this circular orbit we can Taylor expand J{u) around u . Keeping up to first 
order terms in the perturbation of u we get 

J{u) = J{u ) + {u-uo)(^r) ■ ( 10 ) 



du J 



'tin 
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Noting that J(uq) = Uq for the circular orbit, we can now write Eq. (j7|) as 

d 2 u . ( dJ 

+ {u -u ) = (u-u ) 

\ / UO 

If we define x = u — u then the above equation can be written as 

d 2 x 



where ( 2 is defined as 



^+e 2 * = o, (id 



C= = l-f^'l (12) 



\du 



From Eq. ( II ip it is clear that if the orbit of the relativistic particle in a general central 
potential has to be stable then ( 2 > and if the orbit has to be closed then ( must be a 
rational number. 



3.1 A differential equation for the potential V(r) producing sta- 
ble and closed circular orbits 

The rational number ( as predicted, in Eq. (IT2|) . from the stability criterion of closed 
circular orbits in the central force problem is an interesting input in the theory. The 
interesting property about this rational number is that it is a constant and so it does not 
depend on the details of the orbit which one tries to perturb. The reason for the constancy 
of ( is the following. For any circular orbit with radius r$ a specific ( specifies the 
number of undulations of the perturbed orbit. If ( is a rational number then the number 
of undulations of the perturbed orbit will be such that they form a closed geometrical 
structure. Now suppose one takes another circular orbit of radius ro + 5r where 5r <C tq. If 
( has a different value on this orbit then the number of undulations due to a perturbation 
will be different. In the limit 5r — > in a continuous manner the two unperturbed 
circular orbits tends to each other but the number of undulations on the circular orbits 
will not match as ( is not a continuous variable but can only have discrete rational values. 
Consequently the number of cycles of the perturbations will change discontinuously with 
radius and the perturbed orbits cannot be closed at this discontinuity. As we are only 
interested in stable, closed orbits we can conclude that ( must be a constant and not 
change discretely with r. The discussion on the constancy of ( as given above closely 
follows the analysis given in Ref. [T] where the author gives a nice discussion on the role 
of ( in the case of non-relativistic orbits. 

As ( is a constant and must not depend upon the choice of uq or x one can interpret 
Eq. ffl2|) as an independent differential equation by itself, 

whose solutions would give us information about the general form of the central potential 
V(r). Using Eq. (Q we can write the last equation as 
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7 r_2 | ( M = 1 ~ s ; (16) 



which is a non-linear second order differential equation. The right hand side of the above 
equation is a constant which can be written as 

d = L 2 c 2 (l - C 2 ) • (15) 

Eq. (fT4"j) admits multiple solutions for V. The constant E is the total energy of the 
particle. 

It is interesting to note that the differential equation for the potential stemming from 
the stability of closed, circular orbits in the relativistic case does not have a non-relativistic 
analogue. Although the orbit equation Eq. (jSJ) has a proper non-relativistic limit the 
same cannot be said about Eq. (1141) . The reason for such behavior can be seen clearly 
if we rewrite Eq. (fl4"|) in a slightly different way. From the expression of the energy of 
the particle in the central force field 7 can always be written as (E — V)/mc 2 . As the 
total energy is a constant in the present case we must have dj/du = — {\/mc 2 )dV/du. 
Consequently Eq. ([141) can also be written as 

du 2 \du J m 2 c 2 

which gives a differential equation of 7. The equation above obviously does not have 
a well defined non-relativistic limit. The relativistic stability condition produces an ill- 
defined non-relativistic limit due to the fact that in the relativistic case J(u) as given in 
Eq. (j5J) depends upon the velocity of the orbiting particle^. In the non-relativistic J(u) 
was purely a function of the radial coordinate of the orbiting particle. A perturbation 
from the circular orbit in the relativistic case consists of two kinds of perturbations. One 
is related to the change in position of the particle from its previous orbit and the other is 
the change in velocity from the velocity it had previously on the circular orbit. In the non- 
relativistic case only a radial perturbation from the circular orbit fixes the shape of the 
stability condition. As because the stability condition of the orbit depends upon velocity 
of the relativistic particle and the corresponding non-relativistic stability condition does 
not depend upon the velocity of the particle, the non-relativistic limit of Eq. ( 1T4"|) or 
Eq. ([IB"]) is not well defined. 

In the case of non-relativistic motion we know that the inverse square law potential 
and the simple harmonic potential has the capability to produce stable, closed circular 
orbits. In the present case to get the forms of the potentials which can produce stable, 
closed orbits we have to solve Eq. ( I14p . As it is a non-trivial equation we will first try to 
see whether the the potentials which produced stable, closed orbits in the non-relativistic 
regime still satisfy Eq. f JT4|) . Let us try to see whether any power law solution of the form 



V(u) = -au T , (17) 

where a > satisfies Eq. f[T4"j) . In the above equation a and r are constants. If we 
substitute the above form of the potential in Eq. (fi~4"]) we get 

u^-V {a 2 T(r - 1) + a 2 r 2 } - w T ~ 2 {Ear(r - 1)} = d . 

This directly shows that the above relation can be valid for any u only if r = 1, when 
a 2 = d or 

L — a (18) 



1 The V — E in J{u) is proportional to 7 which depended upon the velocity of the particle. 
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on using Eq. ( TT5|) . In this case we see that choosing r = 1 in Eq. ( ITT]) we get the Coulomb 
or Newtonian potential. The last equation shows that for stable, circular orbits the 
particles angular momentum must satisfy some condition. Eq. (ITS]) implies that ( 2 < 1, 
and as ( 2 > for a stable orbit, we have 



The above equation gives 



o < c < i 



r a 

L > - 



(19) 



(20) 



giving a lower bound on the angular momentum of the orbiting particle. This lower 
bound of the orbiting particle was previously obtained in a different way by T. H. Boyer 
in Ref. [2J. It must be noted here that except r = 1 no other values of r are allowed in 
the potential which can produce stable circular orbits of relativistic particles . In non- 
relativistic mechanics we do also have the harmonic-oscillator potential corresponding to 
r = —2 and a < in Eq. f JTTj) . but interestingly relativistic effects forbid this value of r. 



3.2 The general solution of the differential equation for the po- 
tential and the nature of orbits 

We can find out the general form of the force which can produce stable, circular relativistic 
orbits. Noticing that the left hand side of Eq. ( FH|) can also be written as 

du 2 



(V-Ef 



it can be easily shown that 

V(r) 



E 




d b + 



+ a . 



(21) 



satisfies Eq. ( fl4l) where d is as given in Eq. ( fl5l) and b and a are two other dimensional, 
integration constants. For an attractive force b > and d > but a can have any sign. 
If we assume that as r — > oo, V{r) — > then we get a relation between the constants d, b 
and a as 



E = Vdb 2 + a . 
From Eq. (j2T|) we get the force acting on the particle, 

F = -W(r) , 



(22) 



as 



d(b+l) +a 



(23) 



From the form of the force and Eq. (1221) we immediately see that if a = we have 
b = E/\fd and we get back the Newtonian or the Coulombic potential. From the form of 
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the potential as written in Eq. ( ITT]) we can furthermore identify a = \fd and consequently 
when a = we have b = E/a. If a ^ then the form of the force is non-trivial. The 
form of the force as given in Eq. ( 12 3 p cannot be reduced to the harmonic oscillator force 
in any limits of the constants. This shows that special relativistic effects do not allow 
stable circular orbits in the presence of a force which is proportional to the negative of 
the displacement vector. Although the force expression in Eq. ( )23|) is mathematically 
interesting but in physics we do not encounter such a force, except the a = limit. 
From the expression of V — E as given in Eq. ( 12~T]) we get 



(V-E)dV d . 



yielding 



d 2 u 
~dff< 



+ u 



d 



L 2 6 



■{u + b) 



(24) 



(25) 



which gives the orbit equation of the relativistic particle which is acted on by a force given 
by Eq. ( 1231) . As in general d = L 2 c 2 (l — ( 2 ) for a stable closed orbit where ( must be a 
rational number, we get 



I * (ce) + KI_C!) 

r R 



where 



R = Lc( 



b 2 L 2 c 2 (l-C) , 

— ha- me 

C 



-1/2 



(26) 



(27) 



The equation of the orbit in Eq. ([26]) shows that in the most general case we will have 
precession of the orbits dictated by the condition (2tt + 56) £ = 2tt, which predicts that 
the orbit precesses by an angle 



56 



2tt(1 - Q 

c 



(28) 



per orbit. 



3.3 The case of large perturbations 

Till now we have utilized first order perturbation from a circular orbit as described in 
Eq. ([TO]) in the beginning of this section. To include higher order perturbations from a 
circular orbit we require more terms in the Taylor series expansion of J(u) in Eq. ( TTTJT) . 
The second order effects will come from terms proportional to (d 2 J/du 2 ) U0 . If the general 
form of the potential V(r) satisfies Eq. ( |2~T]) then it is immediately clear from Eq. (124"]) 
that all derivatives of J(u), except the first, vanishes. Consequently in the relativistic 
case it is impossible to restrict the constants (, a and b by higher order perturbation 
terms to the circular orbit. For higher order perturbations from circular orbits the form 
the potential as given in Eq. (]2"T]) remains the same. 
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4 Connection of the present work with some related 
works 



One of the findings of the present article is related to the absence of stable, circular orbits 
for relativistic particles in presence of a harmonic oscillator potential. The trajectories of 
relativistic particles in a three dimensional harmonic oscillator potential has been studied 
previously by L. Homorodean in Ref. [BJ. The method followed in the referred work 
is completely equivalent to the one followed in the present work. It is interesting to 
note that in Homorodean's analysis the general shape of the orbit in the relativistic case 
is not an ellipse, or a circle, but a rosette shaped curve. In presence of the oscillator 
potential the angular momentum of the orbiting particle with a specific energy has an 
upper bound. The trajectory of the relativistic particle can only be a circle when it has 
the highest angular momentum for a fixed energy. In Ref. [8] the author does not give 
any information about the stability of the orbits. In the non-relativistic limit the orbit 
of the particle can be circular. In the light of the findings in Ref. [8] of Homorodean the 
prediction of the absence of a stable, circular orbit in the oscillator potential is a sensible 
result. 

Bertrand's theorem in non-relativistic classical mechanics has inspired some authors to 
propose a space-time (a metric to be precise) where any bounded trajectory of a particle is 
periodic in nature. This kind of a space-time is named as Bertrand space-time. The works 
of Perlick, Ballesteros, Enciso, Herranz and Ragnisco, in Refs. [9j [10], try to generalize 
the results of the classical Bertrand's theorem on a flat 3-space to a curved 3-manifold. 
In Ref. [9] the author found that a specific form of a space-time which is asymptotically 
flat can support Keplerian orbits. The asymptotically flat Bertrand space cannot support 
closed trajectories expected in an oscillator type of potential. One of the findings of 
the present article predicts that even in flat space relativistic effects forbid closed, stable 
trajectories of particles in presence of an oscillator potential. 

5 Conclusion 

The outline of the article is based on the well known Bertrand's theorem on central 
potentials and orbits of particles as described in most of the classical mechanics books. 
Like the non-relativistic case the relativistic particle's orbit around a potential source 
takes place in a plane where the angular momentum and presumably the energy of the 
orbiting particle remains constant. The main difference between the non-relativistic orbits 
and relativistic orbits crops up in the orbit equation itself. Unlike the non-relativistic 
case in the relativistic case the orbit equation depends upon the total energy of the 
particle. The main aim of the article was to find out possible forms of central potentials 
which can produce stable circular orbits for relativistic particles. The stability condition 
for the orbits can be transformed to a non-linear differential equation for the central 
potential. It is seen that one of the solutions of the non-linear differential equation for 
the central potential is just the normal Coulomb potential. But relativity affects the 
properties of the orbits by curtailing the angular momentum values beyond a certain 
limit. Except the Coulomb potential solution we find that the stability equation has 
another general mathematically interesting solution which is unlike any potential which 
we use in conventional physics. In a specific limit the general solution reproduces the 
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Newtonian or Coulomb form. In the relativistic version of the central force problem we 
lack some restrictions on the potential which can produce stable circular orbits. In the non- 
relativistic version minimizing the effective potential one can figure out the radius of the 
circular orbits and higher order perturbation corrections to the stability condition of the 
orbits could be used for unravelling the exact nature of the potential. In the relativistic 
version none of those restrictions remain and consequently the general solution of the 
potential contains some constants whose values cannot be analytically calculated. 

An important fact which comes out from the article is about the non-existence of the 
harmonic oscillator potential as a solution of the stability equation. In non-relativistic 
treatment of the Bertrand's theorem it is well known that only two kinds of potentials 
can produce stable circular orbits, one is of the Coulomb type and the other is of the 
harmonic oscillator type. The Coulomb form of the potential passes the stability test for 
circular orbits but the harmonic type does not. 

The article was focussed on some mathematical properties of relativistic particle orbits 
in a central potential. Before we finally conclude it is pertinent to say some thing on the 
practical side of the relativistic central force problem. Atomic physics always remains a 
store house of exciting phenomena and one of the places where one may like to apply the 
tools of relativistic central force problems lies inside the atom. This fact was discussed 
in Ref. [2J. People have studied about the Schrodinger equation and the Dirac equation 
in presence of the Coulomb potential. It can be quite interesting to study the analogous 
problems using the potential presented in this article instead of the Coulomb potential. 
This attempt can seriously shed some light on the physics of the atoms. 

More over as there exists some work on the general relativistic generalization of 
Bertrand's theorem one may expect that in the simplest of the situations, where space- 
time remains flat, the results of the present work can be applied for the orbits of very fast 
moving bodies interacting via Newtonian gravity with a massive source. 
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